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Abstract. We remove a self- interaction from the Galitskii-Feynman T-matrix 
approximation. This correction has no effect in the normal state but makes the 
theory applicable to the superconducting state. It is shown that identical theory 
is obtained by removing nonphysical repeated collisions in the spirit of the Fadeev- 
Watson-Lovelace multiple scattering expansion. Our correction does not violate the 
two-particle symmetry of the T-matrix, therefore the present theory is conserving in 
the Baym-Kadanoff sense. The theory is developed for retarded interactions leading 
to the Eliashberg theory in the approximation of a single pairing channel. 
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1. Introduction 

The family of superconducting materials collects so diverse systems as conventional 
metals and metallic alloys [T], high-T c ceramics [2], fullerenes [3], organic 
superconductors [4], doped diamond J5], heavy fermion metals [6], 3 He [Tj, symmetric 
nuclear matter in nuclei [El EH [TO] and very non-symmetric nuclear matter in neutron 
stars [TT] . Fermi gases [12], not talking about hypothetical condensates like the 
color superconductivity of quarks [13J. It is rather surprising how many features of 
these systems have been successfully explained within the mean-field Bardeen-Cooper- 
Schrieffer (BCS) theory and its Green function extension due to Eliashberg. On the 
other hand, there is a growing list of experimental facts which require to employ more 
elaborate theories. 

Recent theoretical approaches to superconductivity range from trial wave functions 
of Gutzwiller type [H], over improved Eliashberg theories [To] 1TB], renormalization group 
approaches within path integrals [UJ , exact diagonalizations and quantum Monte Carlo 
studies [18j on simple models having small size or infinite dimensions [19] . to the two- 
particle T-matrix approximation [20l ETJ |22} [23] . 

The majority of studies focus on high-order correlations leaving aside details of the 
pairing interaction. Among the most popular approaches, there are the Hubbard model 
and the Fermi gas model with point interaction. Much less attention is payed to the 
retardation and finite range of interaction which are studied only by Eliashberg-type 
theories. For many of systems listed above, however, the retarded nature or finite range 
of interaction are among the most important ingredients of a theoretical model if one 
wants to achieve a quantitative agreement. 

In this paper we derive a T-matrix approximation which combines advantages of the 
Galitskii-Feynman (GF) and the Kadanoff-Martin (KM) approximation. It applies to 
Fermi systems interacting via boson fields and provides the superconducting gap which 
depends on the energy reflecting the retardation of the interaction line. 

2. Self-interaction in the T-matrix approach 

The GF and the KM approximation are compared in a diagrammatic representation 
in Fig. [1] Both approximations are based on the two-particle T-matrix in the ladder 
approximation. As one can see, the KM approximation is nothing more but a simplified 
version of the GF approximation neglecting the exchange and possessing a bare line in 
the closed loop of selfenergy. Nevertheless, these two approximations are very different 
with very disjunct fields of applications. 

The GF approximation is used in nuclear physics for both equilibrium [2H [251 EE] 
and non-equilibrium [271 [28] problems, in the theory of moderately dense gases [29J and 
liquid 3 He [30], and in studies of electron-electron correlations in molecules and solids 
[31~1 1321 [33] |3"H 135] . The KM approximation was never adopted to these problems by 
many reasons. The most important is that the conservation laws are guaranteed only 
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Figure 1. T-matrix approximations in diagrams. Both approximations have a 
selfenergy constructed from the two-particle T-matrix. The interaction carried by 
boson propagators (wavy lines) is included in the ladder approximation. The GF 
approximation has all internal lines of the selfenergy dressed (thick arrows). In the 
KM approximation the loop is made of bare propagators (thin arrows). The GF 
approximation includes the exchange channel in which the lower line matches the 
upper one so that the diagram has no loop. The exchange channel contributes only if 
both particles have parallel spins. 



if the T-matrix is symmetric with respect to the interchange of the upper and lower 
line [36J. The non-symmetry of the KM approximation is viewed as unjustified and 
unacceptable. 

The KM approximation is used exclusively in the theory of superconductivity 
[211 [221 EZl [3HI [391 SOI SU S21 H3l It describes the superconducting gap on 

the level of the mean-field theory and covers the lowest order fluctuations. The GF 
approximation cannot be employed for superconductors in spite of its superiority in 
other fields. Although it becomes unstable at the critical temperature [45J and the T- 
matrix diverges there, the GF selfenergy constructed from the T-matrix fails to describe 
the superconducting gap [16]. 

The paradox that the worse approximation (KM) works well while the better one 
(GF) fails was first noticed by Prange [47J and Wilde [46J prior to the work of Kadanoff 
and Martin. The Prange paradox is not a common knowledge and some authors, see e.g. 
|48j, report a superconducting gap obtained within the GF approximation. By closer 
inspection one finds that they simplify formulas using the bare propagator to close the 
loop. This additional simplification turns the GF into the KM approximation. 

Apparently, the GF approximation includes diagrams which block a formation of the 
gap, while the KM approximation is free of them. One of such contributions is sketched 
in figure [2j The diagram shows two kinds of the nonphysical self-interaction; the simple 
self-interaction in the Hartree approximation and the indirect self-interaction which is 
enhanced by the resonant collisions when the superconducting condensate appears. The 
self-interaction via the Hartree field becomes the important mistake when a particle is 
bounded in the real space. The indirect self-interaction becomes crucial when particles 
are bounded in the momentum space as it is the case of Cooper pairs. 

In Feynman's diagrammatic expansion the self-interaction is canceled by 
corresponding exchange diagrams. Such compensating exchange diagrams are beyond 
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Figure 2. a) Schematic picture of four interacting particles of initial momenta (and 
spins) k, p, q and m. Due to the Pauli principle all particles are in different states, 
e.g. m ^ k and q ^ k. b) Corresponding Feynman diagrams for the Green function 
of the particle k. The other lines are closed into loops and momenta p, q and m are 
summed over with no restriction. For m = k the diagram yields the self-interaction, 
the artifact appearing already in the Hartrcc approximation. For q = k the nonphysical 
self-interaction is more complex. In the superconducting state collisions are enhanced 
for p — —k. Similarly, the enhanced processes couple p with q = —p, which leads 
to q = —p = k. This indirect resonant self- interaction blocks the formation of the 
gap. The g-loop is included in the GF approximation while it is absent in the KM 
approximation. 



the T-matrix approximation and it is not clear how many and how complicated diagrams 
ought to be included. Such cancellation can be demonstrated on the simpler self- 
interaction included in the Hartree field. This self-interaction is canceled by the Fock 
selfenergy - one step more complicated approximation. Apparently, the original simple 
problem is cured only on cost of an increase of complexity of theory. 

One can eliminate the self-interaction restricting sums in the Hartree term in 
specified real space representation. Such correction is not fully equivalent to inclusion 
of the Fock term, but it removes gross mistakes, e.g. the potential binding electrons to 
neutral atom has Coulomb long range asymptotics while in the Hartree approximation 
it is fully screened. Similarly, one can eliminate the self-interaction restricting sums in 
the T-matrix approximation. With respect to Copper pairs, this restriction has to be 
done in the momentum space. This approach we follow here. 

3. Restricted T-matrix approximation 

When a pair of particles interact in stationary and homogeneous systems, its energy and 
momentum Q = (u, Q) conserve. Dressing of a particle of four-momentum k is given 
by the selfenergy E^(k) which is a sum over interacting pairs, H^(k) = J2q °Qt(^)- Here 
cjQ^(fc) is a single Q contribution and J2q ■ ■ ■ = XL Xq • • • denotes sums over bosonic 
Matsubara frequencies and discrete momenta in the quantization volume Q. 

In the spirit of the multiple scattering expansion [19] we introduce a reduced 
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propagator 

GqM = G°(p) + G°(p) £ o-Q^p) G m (p). (1) 

Q'+Q 

The reduced propagator G@i(p) of the particle of four-momentum p and spin |, briefly 
(P'i 4)-particle, is dressed by all binary interactions except for the one in which the sum 
four-momentum of it and its interaction partner equals Q. By this we exclude the 
interaction of the (Q — k; |)-particle with the (k; t)-particle. 

The contribution of a collision of pair four-momentum Q to the selfenergy 

*Qt(*0 = ^r n (k,Q-k;k,Q-k)G m (Q-k) (2) 

we will call a Q-channel of selfenergy for brevity. The dressed propagator G is given by 
the Dyson equation with the unrestricted sum over channels 

G f (k) = G°(k) + G°(k) £ a m (k) G t (fc). (3) 

Q 

To avoid the self-interaction in Eq. (j2J) we use the reduced propagator Gq±(Q — k) 
to close the loop. To guarantee conservation laws, all propagators in the loop have to 
be in the same approximation. Accordingly, the T-matrix has to be constructed as 

T n (k,Q-k;p,Q-p) = D(k,Q-k;p,Q-p) 

_^^D(k,Q-k;k\Q-k')G t (k')G m (Q-k')T n (k',Q-k'-,p,Q-p), 

k' 

(4) 

where D is a bosonic interaction line with interaction vertices included. The prime 
on sum denotes that it runs over fermionic Matsubara frequencies. We follow the sign 
convention of Ref. [50] Sec. 14.2. with D going to the interaction potential in the 
non-retarded limit. 

The set of equations (jX|)- is closed. For the sake of simplicity, we write only 
equations for a selected spin orientation, the complementary equations are obtained by 
flipping all spins, f-< — k|~ These equations do not have diagrammatic representation 
in the strict sense. This is because the Feynman diagrammatic rules are based on 
unrestricted sums over states. 

4. Asymptotic properties 

The theory is complete. Now we discuss its properties. 

In the normal state, the single-channel contribution vanishes in the thermodynamic 
limit, ctq^ oc 1/Q — >■ 0. The reduced propagator then equals to the dressed propagator, 
Gq^ — > Crf, and the present theory is identical to the GF approximation. 

In the well developed superconducting state a single channel becomes singular. In 
the absence of currents and in equilibrium, it is the zero energy and zero momentum 
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channel Q = (0,0). Its T-matrix is separable [5U US] arid diverges being proportional 
to the volume 

T n (k, -k;p } -p) = -JL.A(k)A(p). (5) 

It is advantageous to split the singular part off the remaining regular terms 

S t (A;) = -A(k)G H (-k)A(k) + Epf(k). (6) 

The reminder 

E*t(*) = E *«t(*) (7) 

covers normal processes in the background of superconductivity. From equations §B) 
and (j5D one can see that A does not enter the reduced propagator 

G m (k) = G°(k) + G°(k) E m (k)G m (k). (8) 

With reservations one can say that G^ is a propagator of the normal metal. 

Using the selfenergy (EJ), the dressed Green function (j3J) can be expressed via the 
reduced propagator 

Gf(k) = G m {k) - G m {k) A{k)G n {-k)A{k) G t (k). (9) 

This equation shows that A(k) equals to the energy and momentum dependent 
anomalous selfenergy giving the superconducting gap in the Eliashberg theory [1] . The 
anomalous selfenergy itself follows from the equation for the T-matrix (j4]) and the 
separability (jSJ) 

A(fc) = ~^T! D ( k ^ k',^')G t (k')G H (-k')A(k'). (10) 

k' 

Deriving (TT0|) we have used that D/Q — > in the thermodynamic limit. 

Apparently in the well developed superconducting state the present theory 
asymptotically approaches the Eliashberg theory. There are some differences, however. 
In the present theory all processes, normal and pairing, are treated within the same T- 
matrix approximation. In the Eliashberg theory the normal processes are in the Migdal 
(or Born) approximation while the pairing covered by equations for anomalous functions 
is described by the approximation corresponding to the T-matrix. 

Detail discussion of a region near the critical temperature is beyond the scope 
of this paper. We merely state that the single- channel approximation ( TlOT) can be 
simply modified to non-zero pair momentum giving again the familiar equation for 
the anomalous functions. It implies that phase fluctuations of the gap function are 
the same as those following from the time dependent Ginzburg-Landau theory. There 
are, however, differences in excitation spectrum of noncondensed bound pairs. The 
mean-field theories of Eliashberg type ignore these excitations. The KM approximation 
overestimates their density [52] having quadratic dispersion relation for bound pairs. 
In a next paper it will be shown that the present theory yields an energy gap between 
condensed and noncondensed pairs. 
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5. Two-particle symmetry 

Now we show that the present theory is conserving in the Baym-Kadanoff sense, i.e., 
that it satisfies conditions (A) and (B) from Ref. [36J. The condition (A) links the 
single-particle Green function G with the two particle function Q. For the present 
approximation it is satisfied if the two-particle function relates to the T-matrix as 

G n (k, Q-k; P ,Q-p) = G^k)G m {Q-k) 5{k - p) 

-G t (k)G m (Q-k) T n (k,Q-k;p,Q-p) G t (p)G m (Q-p). (11) 

The condition (B) demands that the two-particle function is symmetric with respect 
to interchange of the upper and lower lines. This symmetry is satisfied although it is 
not obvious from expression ([5]). 

First we show that the T-matrix fll]) is symmetric with respect to the interchange 
of the upper and lower lines 

T n (k,Q-k;p,Q-p) = T^(Q-k,k;Q-p,p), (12) 

in spite of the selfconsistency restricted only in the upper line. Apparently, symmetry 
( TT2|) is satisfied for the lowest order T w D. 

Now we assume that the T-matrix is symmetric in the n-th order and show that 
the symmetry of (n+l)-th order follows. For a general order obtained by iteration of 
Eq. (jl]) it is sufficient to show that 

G m (k)G±{Q-k) = G t (k)G m (Q-k). (13) 

In the iteration process we assume n-th order in the selfenergy of Green functions GqG. 
From Eq. (|2j) we then have 

<T Qt (k)G m (k) = T n (k,Q-k;k,Q-k)G Qi (Q-k)G m (k) 
= T u (Q-k, k;Q-k, k)G m {Q-k)G m {k) 
= G Qi {Q-k)a Qi {Q-k) : (14) 

where we have used symmetry of the n-th order. Now we are ready to proof ( fTB"l) . 
G m (k)G x (Q-k) = G t (k) (1 - a Qt (k)G m (k)) 

x (i + Gi (Q-k)a qi {Q-k))G Qi (Q-k) 
= G t (k)G m (Q-k) (1 - a QX {Q-k)G m {Q-k)) 

x (i + Gi (Q-k)a Qi (Q-k)). (15) 

In the first step we have used 

G^Q-k) = G m (Q-k) + Gi (Q-k)a Qi (Q-k)G m (Q-k) (16) 

which follows from Eqs. (II]) and ([3]), and a similar relation for GQ^{k). In the second 
step we have substituted from Eq. f }T4"|) . From equation f|T6|) follows that the product 
of the two brackets in Eq. ( 1T5|) equals to unity. The relation ( 1T3|) is proved. Therefore, 
the symmetry ( !T2|) of the T-matrix is proved. 
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From the symmetry of the T-matrix (1121) and relation (I13p follows the symmetry 
of the two-particle Green function 

Qn( k ,Q-k;p,Q-p) = Gi,t(Q-k,k;Q-p,p). (17) 

The present theory thus also satisfies the condition (B) and is conserving in the Baym- 
Kadanoff sense. 

We note that in the KM approximation the two-particle Green function satisfying 
the condition (A) is made of one dressed and one bare line. Such two-particle function 
does not satisfy the condition (B). 



6. Multiple scattering approach 

The nonphysical resonant self-interaction is naturally avoided in the Fadeev- Wat son- 
Lovelace multiple scattering expansion |53j EH ES1 ESI EH Ell 09]. In Ref. [23] 
the method of the multiple scattering theory was implemented to the theory of 
superconductivity following a rather different approach. It was argued that the GF 
approximation includes nonphysical contributions in which two particles interact again 
after they have accomplished a collision. Indeed, the T-matrix sums the interaction 
potential to the infinite order so that the next interaction is possible only with a next 
particle. When one eliminates such repeated collisions from the GF approximation, the 
resulting theory goes to the GF approximation in the normal state while it yields the 
gap in the superconducting state. Here we show that the present approach and the 
approach of Ref. [23] modified to four-momentum yield identical results. 

The nonphysical repeated collisions can be eliminated with the help of Soven's 
concept of effective medium modified for binary collisions. In this spirit one selects a 
Q-channel which will be described explicitly while all other channels are described by 
the effective medium represented by the selfenergy. We thus eliminate the Q-channel 
from the selfenergy of the (fc; ^)-particle. This channel is included explicitly into the full 
propagator with a single collision being allowed, 

G t (fc)=G @t (fc) + G 0t (fc)^r n (fc,g-fc;A ; ,Q-A;)G i (g-fc)G 0t (A;).(18) 

Note that the loop is closed by the full Green function. This is because we have 
eliminated ctq^ which does not enter G± in the loop. 

The scattering equation fflBl defines the selfenergy indirectly. Comparing f fl8|) with 
G-f(fc) = GQ^{k) + GQ^{k)aQ^{k)G^{k) we find that Q-channel selfenergy is given by 

i 77s = ^T n (k,Q-k;k : Q-k)G^Q-k). (19) 

1 - (r Qr (k)G m {k) Q 

This relation complicates a straightforward implementation of the theory obtained from 
Soven's scheme. Fortunately, it can be simplified to the set ([T])-(jl]). 

Using the symmetry (Fl3|) one can readily see that equation (fl8|) is equivalent to the 
approximation introduced in this paper. Alternatively, one can multiply Eq. (fl9l) by the 
denominator of the left hand side and use Eq. (|T3j) to turn relation ( 1191 into Eq. (J2]). 
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Finally, we have to take into account that inside the T-matrix the reduced propagator 
appears for spin "f, i.e., in the lower line. According to ffl2|) this is equivalent to the 
T-matrix with the reduced upper line. The present derivation is thus equivalent to the 
application of Soven's scheme. 

We note that here we identify the channel via energy and momentum Q = (u, Q). 
In Ref. [23J the channel was identified only via momentum Q what applies only to 
non-retarded interactions and leads to slightly different results. In particular, the 
identification of a channel via momentum does not provide a two-particle Green function 
symmetric with respect to interchange of the upper and lower lines. The theory in 
Ref. [23] thus does not satisfy the condition (B) and it is not conserving in the Baym- 
Kadanoff sense. 

7. Effect of self-consistency 

The present theory can be implemented to the same family of problems as the KM 
approximation which was recently successfully used to explain many properties of such 
systems as the high-T c materials [371 EH IZH [39] , organic superconductors [4 J and the 
ultra-cold gasses [H E21 SH S21 SSI EH] • 

We do not expect that all results will be identical. As already mentioned, the 
major difference between the present theory and the KM approximation rests in the 
self-consistency with which one includes normal processes in the loop. While the KM 
approximation covers only the Hartree field, the present approximation keeps all normal 
processes described at the level of the T-matrix. 

In this section we demonstrate the effect of the self-consistency on the phase 
diagram of the Fermi gas. For simplicity we use a non-retarded attractive interaction of 
Yamaguchi type 

D(k,Q-k;p,Q-p) = -- £ ^L+pt^Lp, 2 + 2 2 + ^ -14^(20) 

with a soft cutoff at momentum 7. In figure [3] all quantities are dimensionless: the 
temperature is scaled with the cutoff energy e 7 = h 2r ) 2 /2m, the momentum with cutoff 
momentum, k — > k/7, the density with the cutoff cube, n — > n/7 3 . The used interaction 
strength is A — > X/(8nh 2 /Wy) = 10. 

As one can see in figure El at high densities and temperatures the system is a 
Fermi gas or liquid. At lower densities, the Pauli blocking of internal states is not 
effective and diatomic molecules become stable. Region with and without molecules are 
separated by the dash-dot line. The stable molecules appear as a pole of the T-matrix 
at a negative energy - the bound state energy. The grey region at low temperatures 
has non-zero gap. At high densities it corresponds to the BCS condensation of Cooper 
pairs, at low densities to the Bose-Einstein Condensation (BEC) of molecules. The 
dashed line which separates these two regions is the crossover line given by the zero 
value of the chemical potential. Since the chemical potential is fixed to the bound-state 
energy, it can be also interpreted as a line separating stable molecules from Cooper 
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Figure 3. Phase diagram of the Fermi system with attractive interaction. Compared 
to the present theory (thick lines), the Kadanoff-Martin approximation (thin lines) 
reduces the region of the BCS condensate shifting the superconducting phase transition 
towards lower densities. This suppression follows from unequal masses of particles in 
a Cooper pair - an artifact of the KM approximation. 



pairs. The value of the chemical potential is evaluated from the fixed particle density, 
n= k ^fE' k (G t (k) + G i (k)). 

The present theory (thick lines) and the KM approximation (thin lines) yield similar 
phase diagrams, but the transition lines visibly differ. In the KM approximation the 
region of BCS condensate is slimmer. This is a consequence of unequal masses of spin "I 
and ^ particles in this approximation - one mass is dressed while the other is bare. These 
unequal masses complicate a formation of Bogoliubov-Valutin excitations and reduce an 
energy gain due to the gap and also the critical temperature at higher densities. 

In spite of the destructive effect of unequal masses, the KM approximation predicts 
higher critical temperature at low densities. This follows from the employed Yamaguchi- 
type of attractive potential, which reduces the density of states at the Fermi level. In 
the KM theory, the two-particle density of states is only partly renormalized and thus 
higher than in the fully renormalized present theory. 

8. Conclusions and discussion 

In summary, removing the self-interaction from the Galitskii-Feynman approximation 
we have derived a theory which describes the gap including the gap fluctuations. Its 
structure reminds a renormalized Kadanoff-Martin approximation, however, it has four 
major improvements compared to it. First, in the normal state it recovers the well- 
tested original Galitskii-Feynman approximation. Second, the two-particle propagator 
is symmetric with respect to interchange of both lines. This symmetry implies that the 
theory is conserving in the Baym-Kadanoff sense. Third, the removed self-interaction is 
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a correction of the T-matrix which depends on the momentum of the Cooper pair. Non- 
condensed Cooper pairs thus satisfy a different equation than condensed pairs what 
explains stability of the super-currents with respect to excitations of non-condensed 
pairs. Fourth, the present theory is rooted in the standing many-body approach known 
as the multiple scattering expansion. 

We have discussed only singlet channels covered by 7^.. The triplet interaction 
of equal spin components is described by 7^. As long as triplet pairs do not form a 
condensate and no triplet channel is singular, the triplet T-matrix can be handled in the 
GF approximation including the exchange term. We note that the triplet pairing with 
the exchange included is a formidable problem which has not been treated yet within 
the T-matrix approaches to superconductivity. 
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